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ABSTRACT 

We present a study of the kinematics of the isolated spiral galaxy NGC 864, using 
Ha Fabry-Perot data obtained with the GHaFaS instrument at the William Herschel 
Telescope in La Palma, complemented with images at 3.6 /im, in the R band and in 
Ha filter, and integral field spectroscopic data. The resulting data cubes and velocity 
maps allow the study of the kinematics of the galaxy, including in-depth investigations 
of the rotation curve, velocity moment maps, velocity residual maps, gradient maps 
and position- velocity diagrams. We find asymmetries in the velocity field in the bar 
zone, caused by non-circular motions, probably in response to the potential of the 
bar. We also find a flat-profile bar, in agreement with the strong bar, with the grand 
design spiral pattern, and with the gap between the ends of the bar and the start of 
the spiral arms. We quantify the rate of massive star formation, which is concentrated 
in the two spiral arms. 

Key words: galaxies: kinematics and dynamics - galaxies: spiral - galaxies : individ- 
uals : NGC 864 



1 INTRODUCTION 

Galaxies are the basic building blocks of the Universe and 
their formation and evolution is of great interest in current 
astrophysical research. Their dynamics and morphology are 
* Email: serroz@iac.es the result of both externally driven (e.g. galaxy mergers) 
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and secular evolution (e.g. bar or spiral pattern driven). 
Disentangling the different early evolutionary tracks is com- 
plicated by redshift, distance and dust absorption. Alterna- 
tively one can explore the fossil record of evolution through 
detailed observations of the end product, nearby galaxies. 

One of the most important drivers of this internal sec- 
ular evolution is the flow of gas into the central regions. 
This is the result of angular momentum loss in shocks in- 
duced by non-axisymmetric potentials, as generated by a 
bar, interactions, minor mergers or even by minor deviations 
from axisymmetry such as ovals, lenses, lopsidedness or spi- 
ral arms (e.g., [Schwarz|1984| |Shlosman, Frank &; B egelman 
1989) IShlosman, Begelman fe Frank||1990| |Rix fe Zaritsky| 
T^|Knapen et al.||1995| |Kormendy fe Kennicutt||2004l 



(e.g.,|Marcelin, Boulesteix &; Georgehn|1985||Bonnarel et al. 
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Comeron et al. 2010). Around two thirds of local galaxies 
exhibit at least one bar ( de Vaucouleurs et al.|T99l| Sell- 
wood &: Wilkinson| [1993| IMoles Marquez &; Perezl 1995 1 
Ho, Filippenko fe Sargent |1997[ |Mulchaey fe Reganjp^n 
Hunt fe Malkan||1999| |Knapen, Sh losman fe PeletierPOOol 



Eskridge et al. 2000| |Laine et al. |2002| |Laurikainen, Salo 



ButaM2004l 
fc Jogee|2007] 



Menendez-Delmestre et al.||2007[ |Marinova 
Sheth et al.||2008| |Laurikainen et al.||2009| ), 



so potential drivers of secular evolution are common. Low- 
mass disks acquired their bars later than more massive ones 
(She th et al.|2008 ), implying that the bar-driven secular evo- 
lution is expected to be different depending on each galaxy's 
mass. 

To study the influence of the past evolution of a galaxy 
on the observed morphology, kinematic information is es- 
sential. The study of velocity fields in spirals has mostly 
been done using the 21 cm HI line, primarily because this 
emission can be traced far out, often three or four times 



beyond the visible disk. Rots (1975), van der Hulst (1979), 



Bosma (1981), Gottesman (1982) and others demonstrated 



the power of this technique in deriving the total mass dis- 
tribution of disk galaxies. Since the beam sizes associated 
with 21 cm observations were typically at least as large as 
the disk scale length of the galaxy, little small-scale struc- 
ture in the velocity field was detected in this manner. The 



HI Nearby Galaxy Survey (THINGS; [Walter et al.| |2008) 
presents some of the highest spatial resolution (~ 6 ) 21-cm 
HI observations of nearby galaxies to date, using the Very 
Large Array (VLA) of the National Radio Astronomy Ob- 
servatory (NRAO). 

The best candidate technique to yield two-dimensional 
kinematical maps across a whole galaxy with good arcsec- 
onds spatial resolution is Ha. Due to the cosmic abundance 
of hydrogen, Ha is often the brightest emission line in the 
visible wavelength range. In spiral galaxies, this line traces 
primarily the ionized gas in HII regions around young mas- 
sive stars. Traditional long slit spectra with optical emis- 
sion lines are valuable for deducing the rotation curve of a 
galaxy (e.g., Rubin, Ford & . Thonnard 1980 | |Rubin et al" 
T982l|Rubin et al.|19851|Mathewson Ford &; Buchhorn|19921 
Amram et al.||1992[ |Amram et al.||1994[ or jMathewson fe| 
Ford|1996| ), and the rotation curves can be extended to very 
large radii using Ha (e.g., |Christlein fe Zarit sky 2008). Ha 
observations using Fabry-Perot (FP) instruments have been 
used for some 40 years now (|Tully||1974||Deharveng fe Pel- 
let] [19751 |Dubout et"aLl|1976| or |de Vaucouleurs fe Pence 
1980), and have since been employed to create high signal- 
to-noise ratio (S/N) rotation curves for many spiral galaxies 



1988[ [Pence, Taylor fe Atherton|1990| |Corradi fe CapaccioTi 



1991 JCecil, Wilson fe Tully 19 92) |Sicotte, Carignan fe Du- 
rand||1996[ |Ryder et al.||1998[ |Jimenez- Vicente et aL 1999; 



|Garrido et al.||2002[ |Epinat, Amram k, Marcelin||2008| and 



many more). 

Using the Infrared Array Camera (IRAC; |Fazio et al, 



2004) operating at 3.6 and 4.5 /im, the Spitzer Survey of 
Stellar Structure in Galaxies (S 4 G; Sheth et al.|2010| targets 
over 2300 galaxies. The sample selection criteria is volume-, 
magnitude- and size- limited (<i<40 Mpc, mb<15.5, D25>1 
arcmin). The cornerstone of the survey is the quantitative 
analysis of photometric parameters, enabling a variety of 
studies on secular evolution, outer disk and halo formation, 
galaxy morphology etc. The new images, much deeper than 
traditional ground-based near-IR observations, will allow a 
comprehensive and definitive study of galaxy structure not 
only as a function of stellar mass, but also as a function of 
environment, vital to test cosmological simulations predict- 
ing the mass properties of present day galaxies. 

We have designed an observing programme to obtain 
Ha FP kinematic data sets of over 40 spiral galaxies that 
are included in the S 4 G sample. The observation of NGC 
864 is part of this wider programme using the instrument 
Galaxy Ha Fabry-Perot System (GHaFaS) at the William 
Herschel Telescope (WHT) in La Palma. NGC 864 is used 
to illustrate the data and methods, and discuss the kinds of 
results the main survey will give. 



2 TARGET SELECTION 

The galaxies in our FP programme had to satisfy the follow- 
ing requirements: first of all, the galaxy should fit well in the 
GHaFaS field of view (FOV) of 3.4 x 3.4 arcmin. Therefore, 
galaxies with diameters between 2.7 and 3.4 arcmin were se- 
lected. The declination should be higher than -10° so that 
the altitude in the sky can be enough at the time of the 
observation in La Palma, the Hubble type should be later 
than Sab (T=2; RC3, |de Vaucouleurs et al.|1991 ) to ensure 
Ha emission, and only galaxies with inclinations between 0° 
and 70° were selected. 

These criteria yielded a sample of 108 S 4 G galaxies, 
from which the observed sample was selected by requiring 
a spread in morphological type, galaxy mass (from absolute 
magnitude; ten galaxies in each of four bins), bar presence, 
visibility on the sky and finally, the availability of ancillary 
data (interferometric CO, HI, UV, Spitzer MIR or Herschel). 
The galaxies observed in our first observing run were NGC 
428, NGC 864, NGC 918, NGC 1073 and NGC 7479. This 
paper will concentrate on the galaxy NGC 864. 

NGC 864 is a spiral galaxy, well isolated from nearby, 



similar-sized companions (CIG, |Karachentseva 1973), al- 
though it is a radial-velocity confirmed member of a small 
group (Fouque et al. 1992), including UGC 1670 (v =1601 
km/s; AM = 3.12 mag; ~ 84 arcmin or ^5.7 Mpc away 
from NGC 864) and UGC 1803 (v = 1624 km/s; AM 
— 3.0 mag; ~89 arcmin or ~6.1 Mpc away from NGC 
864); all these values taken from NASA/IPAC Extragalac- 
tic Database (NED). NGC 864 is cla ssified as SAB(rs)c by 
The de Vaucouleurs Atlas of Galaxies ( |Buta, Corwin &; Ode- 
|wahn|[2007| . |Buta et al.| ( [2005] ) found a maximum relative 
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Table 1. Galaxy global properties. (1) The centre position was 
fixed to the position of the optical centre from jLeorT fe Verdes- 
Montenegro (2003). (2) The systemic velocity is listed as the cen- 
tral velocity of the HI spectrum measured at the 20% level, from 
Espada et al. (2005). (3) Distance calculated after applying the 
Virgo, GA and Shapley corrections, with Ho= 73 ± 5 km s _1 
Mpc~\ from NASA/IPAC Extragalactic Database (NED). (4) 
Morphological type from The de Vaucouleurs Atlas of Galaxies, 
( |Buta, Corwin fc Odewahn|[2007| . (5) B magnitude f rom The 
Third Reference Catalogue of Bright Galaxies (RC3; |de Vau^] 
couleurs et al.|1 991). (6) Bar and Spiral Arm grow torques, from 



of the kinematic major axis, from the RC3. (8) Disk inclination 
and PA obtained from the 25.5 mag/arcsec 2 isophote with the 
S 4 G 3.6 jum image, calculations and uncertainties explained and 
to be published in Muhoz-Mateos et al. (2012, in prep). (9) Disk 
inclination and PA from our Ha FP data (Section 4.2.1). (10) 
Luminosity and Star Formation Rate (SFR) of the whole galaxy 
measured from our H a images (Section 4.3). 



Name 



NGC 864 



Comments 



«J2000.0 


2^15 m 27.6 s 


(1) 


<5j2000.0 


+06°00 / 09.1 // 


(1) 


-Usys (kms -1 ) 


1561.6 


(2) 


D (Mpc) 


20.9 ± 1.5 


(3) 


Type 


SAB(rs)c 


(4) 


m B (mag) 


11.62 ± 0.21 


(5) 


Qb 


0.321 


(6) 


Qs 


0.134 


(6) 


i>RC3 


43° 


(7) 


^3.6/i.m 


44.6° 


(8) 




28.2°± 20.4° 


(9) 


PA^cs 


20° 


(7) 


PA3. 


22.0° 


(8) 




25.3 °± 6.7° 


(9) 


L(Ha) (erg s" 1 ) 


(4.0 ±1.6) • 10 41 


(10) 


SFR (M yr- 1 ) 


2.19 ± 0.88 


(10) 



bar torque Qb=0.321, implying that the bar is fairly strong 
in spite of its apparent weakness in blue light. The main 
characteristics of the galaxy are presented in Table [I] 

This galaxy has been observed and analysed by many 
other authors, whose work we summarise here. Ma rtini et al.| 
( 2003 ) studied visible and near-IR Hubble Space Telescope 



(HST) images and color maps of NGC 864, identifying bright 
star formation rings. One of these is a small circumnuclear 
star-forming ring with a radius of less than 100 pc, confirmed 
by the Atlas of Images of NUclear Rings (AINUR; |Comer6n| 
|et al.||2010| ). |Epinat, Amram fe Marcelin| ( |20Q8| ) described 
the analysis of Ha FP cubes of 97 galaxies included in their 
GHASP survey of 207 spiral galaxies, presenting the rotation 
curve of every galaxy, including NGC 864. |Buta & Zhang 
(2009) found three corotation radii in the galaxy: one lies 
within the bar, the next is located slightly beyond where 
the bar ends, and the third is within the inner spiral arms. 



Espada et al. ( 2005 ) provide the most extensive and sig- 
nificant study of NGC 864, via HI observations. They mea- 
sure large-scale asymmetries in the two-dimensional kine- 
matics, finding that the galaxy's profiles are symmetric in 
velocity but asymmetric in intensity. These asymmetries are 
found in the outer parts of the disk, beyond the extent of 
our new data. 



3 OBSERVATIONS AND DATA REDUCTION 

The data for NGC 864 have been taken at different epochs 
using three different instruments of the 4.2-m WHT on the 
Roque de Los Muchachos Observatory in La Palma: the 
GHaFaS Fabry- Perot instrument, the AC AM camera and 
the SAURON integral field spectrograph. These data are 
complemented with a 3.6 /xm Spitzer image from the S 4 G 
survey. 



Buta et al.| ( |2005| >. (7) Disk inclination and position angle (PA) 3 ^ AC AM 



One of the motivations of this paper is to develop a method 
to flux-calibrate GHaFaS Fabry-Perot data cubes. To do 
this, we have obtained narrow-band imaging of the galaxy 
nearly simultaneously with the GHaFaS data. NGC 864 was 
observed on 2010 September 4 with the ACAM instrument, 
permanently mounted at a folded-Cassegrain focus of the 
WHT ( |Benn, Dee fe Ag6cs|2008| . The FOV in imaging mode 
is 8 arcmin, with a pixel size of 0.25 arcsec. 

Two images were taken with ACAM: a 30 s exposure 
image using a Johnson filter R, with central wavelength 6228 
A and FWHM 1322 A; and an Ha image, with an exposure 
time of 300 s, using a narrow filter having a central wave- 
length 6589 A and FWHM 15 A. The [NII]-6548 line falls 
into the narrow filter transmission range (at 6582 A), with 
a similar throughput as the Ha line (which lies in 6597 A). 
However, the contribution of the [Nil] line is not important, 
as in the majority of star forming regions in galaxy disks the 
[Nil] fraction is small or negligible (Jam es et al.|2005 ). 

The basic reduction of these images was performed 
mainly using IRAF. Firstly, bias and flat corrections were 
made. Then, the continuum was subtracted following the 
procedures outlined in |Knapen et al.| ( |2004| and the image 
was photometrically calibrated. To perform the continuum 
subtraction, the i?-band image was scaled by a factor of 
0.0109 and subtracted from the Ha image. The resulting 
images are shown in Figure [I] 



3.2 GHaFaS 

3.2.1 Observations and data reduction 

The galaxy was observed for 4 hours on the same night as 
the ACAM images with GHaFaS ( |Carignan et al.||2008| ), 
situated in the Nasmyth focus of the WHT. GHaFaS is a 
Fabry-Perot imaging interferometer that provides high spec- 
tral resolution Ha images with a seeing-limited spatial reso- 
lution within a 3 '4 x 3 .'4 FOV. The instrument is composed 
of a focal reducer, a filter wheel, a Fabry-Perot etalon, an 
image photon-counting system (IPCS) camera, and a cali- 
bration lamp (neon source). The etalon employed works at 
an interference order of 765 and a Free Spectral Range (FSR) 
of 391.9 km s _1 . Each observational cycle consists of step- 
ping through the etalon for 10 seconds per channel. In the 
observation of NGC 864, 29 cycles of 48 channels/cycle were 
observed in total. The high spatial resolution mode was used, 
achieving a 8.167 km/s velocity resolution with O'.'2/pix scale 
in lKxlK pixel images. 

As a proper de-rotator at the Nasmyth focus of the tele- 
scope is not provided for GHaFaS, the Fabry-Perot cubes 
need to be de-rotated before being further reduced. To 
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Figure 1. i?-band (top) and continuum-subtracted Ha (bottom) 
images of NGC 864, taken with ACAM. 



do this, two fixed points in the field are selected in ev- 
ery image (ideally stars), and used to de-rotate each plane. 
Lamp exposures are used to calibrate the observations in 



wavelength (see Carignan et al. 2008). Due to the nature 
of the Fabry-Perot data, any two-dimensional transforma- 
tion (translation or rotation) of the raw images has to be 
done three-dimensionally (including the third spectral di- 
mension), and thus wavelength calibration and de-rotation 
must be performed simultaneously. The resulting de-rotated 
and wavelength-calibrated data cube is then imported into 
GIPSY (Groningen Image Processing System; van der Hulst 
et al. 1992). The data are placed on an astrometrically cor- 
rect spatial grid by comparing the positions of stars in the 
ACAM images, yielding a data cube of 48 planes of 1001 x 
1001 pixels each. A very bright star nearby, approximately at 
RA 2^15 m 30 s and DEC +05°59 / 49 // , was masked out with 
the GIPSY task BLOT . 

De-rotating blurs the image somewhat, and the FWHM 
after the de-rotation was determined to be 1V56 along the 
RA-axis and l'/4 along the DEC-axis (the worst seeing that 
night was around l'/2, with stages of 0'/8 during the 4 hours 



of observation). The cube was smoothed to different reso- 
lutions with the GIPSY task SMOOTH. We have worked 
at two different spatial resolutions. A first slight smoothing 
was done so that the image loses very little spatial resolution, 
hence taking advantage of the high spatial resolution of the 
instrument (high-resolution cube hereafter, with a FWHM 
of l'/65 x 1V56). Then, a 3" x 3" smoothing was performed in 
order to get continuous spatial information and to highlight 
low-emission regions, facilitating the computation of, for in- 
stance, the rotation curve (hereafter low-resolution cube). 
After the spatial smoothing, the continuum was subtracted 
using CONREM of GIPSY, which estimates the continuum 
level from the line-free frames (the first three channels and 
the last seven in this case). 

To create moment maps, the GIPSY task MOMENTS 
was used. The moment maps start from the smoothed and 
continuum-subtracted cubes. We imposed two conditions: 
the signal has to be higher than 3cr, in at least three ad- 
jacent channels. The resulting moment maps are presented 
in Figure [2] There is a little leakage from the neighbour- 
ing [NIIJ-6548 line, that has to be taken into account when 
calibrating the flux. However, whenever a second peak is 
present, the task MOMENTS does not take it into account 
because it only uses the main peak. 



3.2.2 Flux calibration 

This paper presents the first attempt to flux- calibrate a 
GHaFaS data cube. Previous observations with the instru- 
ment were used for kinematic studies and therefore didn't 
need being flux-calibrated. The continuum-subtracted and 
flux- calibrated ACAM Ha image was used to perform the 
flux calibration of the GHaFaS data. Fluxes from selected 
HII regions in both the GHaFaS cube and the ACAM im- 
age were measured, and then compared. The transmission of 
the filter has to be taken into account, assuming a Gaussian 
profile. As explained before, the ACAM Ha image contains 
a little leakage coming from the [NIIJ-6548 line, and the 
FP also used a filter that allowed light from that [Nil] -line 
to come in. It is thus necessary to correct for these leaks, 
and the correction applied is of 3% following the procedures 
of |James et al. ( 2005 ) . Figure [3] presents the linear fit to 
data points corresponding to a set of HII regions from the 
GHaFaS and ACAM data. 

The uncertainty in the flux calibration is 10% to 20%, 
arising mainly from the use of galaxy distance in the cali- 
bration of the ACAM image, and to a lesser degree from the 
method described above. 



3.3 SAURON 

NGC 864 was also observed with the two-dimensional 
integral-field spectrograph SAURON (Spectrographic Areal 
Unit for Research on Optical Nebulae, |Bacon et al.||2001|), 
and th ese data of NGC 864 were presented in Ga nda et al.| 
( |2006| ). SAURON delivers a spectral resolution of 4.2 A 
FWHM, spatial sampling of 0.96 arcsec and covers the nar- 
row spectral range 4800-5380 A (1.1 A pixel" 1 ). This wave- 
length range includes a number of important stellar absorp- 
tion lines (e.g., H/3, Fe, Mgb) and potential emission lines as 
well (H/3, [OIII], [NI]). The field of view, of 33 x 41 arcsec 2 , 
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Figure 2. a) Intensity map for the high-resolution data. The centre is marked with A. B is the brightest HII region, located at the end 
of the western arm. C and D refer to other regions of massive star formation, found at the ends of the bar, and near the beginning of 
the arms, b) Velocity map for the high-resolution data, c) Dispersion velocity map for the high-resolution data, d) Intensity map for the 
low- resolution data, e) Velocity map for the low-resolution data, f) Dispersion velocity map for the low-resolution data. The centre of 
the galaxy is indicated with a cross, calibrated colours are labelled in the vertical colour bars. 
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Figure 3. Linear fit for the flux calibration. With squares, fluxes 
of some HII regions measured in ACAM Ho; image presented ver- 
sus the fluxes of the same regions measured in the GHaiFaS Ha 
intensity map, whereas the straight line represents the line fitted 
to the data. 



is much smaller than those of GHaFaS or ACAM, but the 
instrument delivers stellar as well as ionized gas kinematic 
information simultaneously. 



The data for this galaxy were taken in 2004 January 20-26. 
Four exposures of 1800 s each were obtained using the low 
spatial resolution mode of SAURON, with the FOV men- 
tioned before. 



Ganda et al. (2006) reduced the data using the ded- 



Ganda et al. (2006) presented the stellar and gas kine- 



matics of 18 nearby late-type galaxies, including NGC 864. 



icated software XSAURON, including bias and dark sub- 
traction, rotation of all the frames, extraction of the spectra 
using the fitted mask model, wavelength calibration, low- 
frequency flat-fielding and cosmic ray removal. Ater that, 
the data reduction consisted of homogenization of the spec- 
tral resolution over the FOV, sky subtraction, and finally 
flux calibration of the spectra, although the data were not 
necessarily collected under photometric conditions. In Fig- 
ure [4] the velocity maps for gas (ft/3) and stars (computed 
using the wavelength region that contains the H/3, Fe5015 
and Mgb absorption lines; see Ganda et al. 2006 for more 
details) are presented. 



3.4 S 4 G image 

In order to derive geometric galaxy parameters and to study 
specific galaxy features, a 3.6 /im image from the S 4 G sur- 
vey has been used. This image was observed in the warm 
mission phase of Spitzer. As the galaxy diameter is lower 
than 3'3 and the Spitzer FOV is 5', the galaxy was mapped 
using a single pointing, following a standard cycling small 
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Figure 4. a) Velocity map of the stars in NGC 864, observed 
with SAURON. b) Velocity map of the gas (H/3) observed with 
SAURON. The kinematic minor axis (PA=20) is plotted with a 
dotted line, whereas the straight line shows the derived systemic 
velocity (see Section 4.2.1). Center is indicated with a cross. 



dither pattern with 4 exposures of 30 seconds each in two 
separate Astronomical Observation Requests (AORs). The 
image has a pixel scale of ; /75, with an angular resolution 
of l'/7. The reduction pipeline is described in detail in |Sheth| 
et al.l (12010). 



4 RESULTS 
4.1 Morphology 

NGC 864 is a late-type barred spiral galaxy, with strong 
emission in the centre. In blue light, the type is SAB(rs)c 
(The de Vaucouleurs Atlas of Galaxies; Buta, Corwin fe| 
|Odewahn|2007t , but the mid-IR type is more like SAB(rs)b 
or be. The bar does look stronger in the S 4 G image (Figure 
|5j), but not as strong as in type SB. The bar is foreshort- 
ened, which undoubtedly contributes to its weak projected 
appearance. 

The bar and some dust lanes can be distinguished in the 
ACAM R-b&nd image (Figure [I]). The bar presents bright 
regions of massive star formation at the ends, where two 




Figure 5. Top: 3.6 /im S 4 G image. Bottom: the same with Ha 
contours at 1.1 ■ 10 36 erg/s, taken from the ACAM image. 



asymmetric arms arise. The arms do not, as is usually the 
case in galaxies, start at the end of the bar, but in fact 
overshoot that position, leading to a mismatch between the 
start of the spiral arms and the ends of the bar. We measure 
from the S 4 G image that the eastern arm starts 0.89 ± 0.07 
kpc away from the east end of the bar, while the western 
arm arises 0.96 =b 0.07 kpc from the west end of the bar. 
This gap between the end of the bar and the start of the 
spiral arms may be the space between the inner resonances 
and corrotation, as spiral arms may be excited only beyond 
corrotation (like |Elmegreen &; Elmegreen|1 985 found in some 
flat-profile bar spirals). 

The arms make less than a half-turn before disappear- 
ing. The outer spiral pattern is more flocculent than a 
two-armed spiral, but not as flocculent as, e.g., NGC 2841 
(Block, Elmegreen &; Wainscoat|1996|. This was also noticed 



by Elmegreen &; Elmegreen| ( |1987| , as they classified NGC 
864 to be a grand design galaxy with Arm class 5, that is, 
two symmetric, short arms in the inner regions and irregular 
outer arms. 

The skewed appearance of the bar in NGC 864 seems 
to be a result of projection effects. It is possible to distin- 
guish both the outer flat part of the bar (i.e. the part of the 
bar coplanar with the disk) and the inner vertically thick- 
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ened part of the bar (usually referred to as a box/peanut 
'bulge'). Athanassoula & Beaton (2006) presented N-body 



simulations of barred galaxies explaining how projection ef- 
fects can lead to the same configuration in the case of M31. 
In NGC 864, the flat outer part of the bar has a position 
angle of about 90-100 degrees (counter-clockwise from N), 
whereas the inner part of the bar has a position angle of 
around 70-80 degrees. 

Alternatively, the skewed appearance of the bar might 
be real. If we inspect the bar starting from the galaxy cen- 
ter, it looks straight initially and then curls a little like 
an open spiral, although the overshoot of the main arms 
could contribute to this. In addition, the inner section looks 
like a twisted oval, and not much like a box with a subtle 
X-pattern. In the de Vaucouleurs Atlas of Galaxies ( |Buta, 



Corwin & Odewahn 2007), clearer cases of low inclination 



examples of this are present. 

We have performed a bulge/disk/bar 2D decomposi- 



tion of the S 4 G 3.6/im image of NGC 864, using budda dde 
Souza, Gadotti fe dos Anjos||2004[ |Gadotti|[200^ ), in order 



to check whether the light profile along the bar is exponen- 
tial or rather flat. [Elmegreen & Elmegreen] \ 1 985| found that 
bars in galaxies with morphological types later than Sc show 
exponential light profiles, whereas bars in galaxies with ear- 
lier types show flat profiles. The RC3 classification of NGC 
864 is Sc, and thus it is clearly a borderline case. The fit is 
complicated by the projection effects mentioned above and 
the rather unusual start position of the spiral arms. We find 
that the Sersic index of the bar is 0.24, meaning that the bar 
profile is flat, and even somewhat brightening at the ends, 
although this is likely an effect of an enhanced luminosity 
at the ends due to a pair of bar ansae and the beginning 
of the spiral arms. Nevertheless, the bar profile is certainly 
not exponential. To verify the robustness of these results, 
several fits were performed with varying initial parameters. 
Although the results of each of these fits show variations for 
a few parameters, the bar Sersic index is always lower than 
0.5, confirming a flat light profile, with high confidence. 



Elmegreen et al. (2011) present a study of spiral arm 
properties of 46 galaxies in the S 4 G sample, which in ef- 
fect forms a modern-day extension of the pioneering work of 



Elmegreen & Elmegreen (1985). In their sample, all of the 



flocculent barred galaxies have exponential bars, while all of 
the grand design barred galaxies have flat bars. NGC 864 is 
a grand design galaxy with a flat profile and thus fits into 
this picture. 

Athanassoula & Misiriotis (2002} showed that flat pro- 



files are found in strong bars and exponential ones in weaker 
bars. This fits in with the results as presented by [Combes 
&; Elmegreen (1993), who found from simulations that flat 



and exponential bars are linked to the shape of the rotation 
curve and halo concentration, finding that flat bars tend to 
occur in galaxies with sharply rising inner rotation curves 
(as we find in NGC 864 - see Section 4.2), while exponential 
bars tend to occur in galaxies with slowly rising rotation 
curves. Athanassoula (2003) showed that this is due to the 



fact that sharply rising rotation curves are linked to a strong 
spheroidal component in the inner parts, and a considerable 
amount of spheroidal mass around the main resonances. This 
leads to significant angular momentum exchange within the 
galaxy, which in turn can lead to strong bars. This was fur- 



the two bar types and compared them to observations made 
by |Gadotti et al.| ( |2007) . 

NGC 864 is a clear case where a flat-profile bar is found 
in a two-armed spiral with a strong bar. The fact that there 
is a gap between the ends of the bar and the beginning 
of the spiral arms is also consistent with flat-profile bars 
( Elmegre en &; Elmegreen|1985| . The bar does indeed appear 
to be of the ansae type, which is very unusual for a late- type 



galaxy ( Martinez- Valpuesta, Knapen &; Buta 2007). This 
means that NGC 864 is an example of an early-type bar in 
a late- type galaxy. Its spiral structure, however, is not very 
grand design even in the mid-IR. Also, it is not an ansae 
type in the typical manner since it looks like an open spiral. 

The moment maps (left plots in Figure [5| and the 
ACAM Ha image (Figure [T| highlight the regions with 
the most Ha emission. We note two main emitting regions 
marked in Figure [2] panel a: the centre (marked with A), 
and one the brightest HII region located at the end of the 
western arm (letter B). The regions where the arms start 
and the bar ends have also been marked with C (east) and 
D (west). Such morphology is common in barred galaxies. 



4.2 Kinematics 

We analysed moment maps at high and low resolution from 
the FP GHaFaS observations, and gaseous and stellar mo- 
ment maps are from the SAURON observations. 

The patchiness of the Ha emission maps doesn't al- 
low us to analyse the maps without smoothing, so for a 
qualitative study of the overall velocity fields, the low- 
resolution (3" x 3") smoothed map is used. However, the 
high-resolution map allows us to study the kinematics of 
the faintest regions. 



4-2.1 Rotation curve 

The rotation curve of the galaxy can be derived from the 
velocity map, which contains the projected velocity along 
the line-of-sight (vi os ) at each position of the galaxy. Only 
LOS velocities are observable so it is necessary to correct for 
projection effects. 

The standard method for deprojecting galaxies from 
kinematic data is to assume that the galaxy is made up of a 
number of annuli, the well known tilted-ring model ( |Bege-| 
man|198"9t , where each ring can be defined with the parame- 
ters: inclination (i), position angle (PA), centre (xq and yo) 
and the systemic velocity (v sys ). 

Rotation curves have been obtained with ROTCUR 



procedure in GIPSY, following the Begeman (1989) method 



Most initial entry values were taken from the literature, as 
in Table ^ but the initial rotation velocity values have been 
taken from the rotation curves in Espada et al. (2005). Set- 



ther discussed by Athanassoula et al. (2009), who modeled 



ting the expansion velocity to zero always, the following 
steps were carried out to obtain some galaxy parameters 
(galaxy centre, PA, i) and the final rotation curve. (1) All 
parameters were left free. The resulting rotation curve was 
not correct and some parameters had to be fixed. (2) To ob- 
tain the position of the galaxy centre, the PA, inclination, 
systemic and rotation velocities were fixed, and the centre 
values were free. (3) The centre position, PA, inclination, 
and rotation velocities were fixed in order to fit the systemic 
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Figure 6. Rotation curves for NGC 864. The black line corre- 
sponds to the rotation curve derived from the GHaFaS data, using 
the low-resolution velocity map. The dark blue plus signs show 
the rotation curve derived from the SAURON H/3 velocity map, 
whereas the red diamods show the rotation curve derived from the 
SAURON stellar velocity map. The green triangles correspond to 
the rotation curve from HI observations taken from Espada et al. 
(2005). The pink squares and light blue crosses represent the ro- 
tation curves for the approaching and receding halves of NGC 
864 derived in Epi nat, Amram fe Marcelin| ( |2008| >. Notice that 
the FOV of SAURON covers only the inner part of the galaxy. 



velocity. (4) The PA and inclination were then left free with 
all other parameters fixed. (5) Finally, the rotation curve 
was obtained by fixing all parameters except for the rota- 
tion velocities. The final values for PA and i are presented in 
Table 1, whereas the resulting rotation curves can be found 
in Figure [6] Notice that the resulting rotation curves take 
into account the deprojected velocity using the computa- 
tion of the different parameters. The same procedure has 
been followed to obtain the rotation curves derived from the 
SAURON stellar and gaseous velocity maps (Figure [4]), and 
the results are also presented in Figure [5] 

Prior to discussing the physical implications of the fea- 
tures in the observed rotation curve, we briefly discuss the 
rotation curves of NGC 864 found in the literature (from HI 
observations, Espada et al. 2005, reproduced in Figure[6]un- 
til 120 arcsec; and from Ha; observations, |Epinat, A mram & 
|Marcelin|2008| as compared to the rotation curves obtained 
using the GHaFaS FP data and the SAURON data. It is ob- 
vious that the information provided from each rotation curve 
is different, as each wavelength offers different information at 
different radial angles. On the one hand, the rotation curve 
derived from HI observations suffers from beam smearing 
due to the low spatial resolution, although it reaches the 
outer parts of the galaxy disk. On the other hand, rotation 
curves obtained from Ha observations are based on detailed 
but also more patchy emission, which can make the analysis 
more difficult. 



The complete HI rotation curve (Figure 7 in Espada 
et al.|2005 ) shows the behaviour that is typical of a galactic 



disk: solid-body rotation near the centre, and almost con- 
stant velocity in the outer parts. At approximately 100 arc- 
sec from the centre (~ 10 kpc), the curve reaches a maximum 
of 130 km s _1 , and continues more or less at that velocity 
until the observed disk ends (~ 40.5 kpc). 



The rotation curves derived from SAURON data show 
the first 20 arcseconds from the center, which presents essen- 
tially solid-body rotation. However, the rotation curves de- 
rived from Ha observations extend as far as the optical disk, 
some 100 arcseconds from the center (~ 10 kpc). The cen- 
tral parts show a quick rise in velocity, which is followed by 
an essentially flat rotation curve for R > 25". Although the 
receding and approaching halves have been analysed sepa- 
rately in the Ha rotation curves found in the literature (pink 
squares and blue cr osses in Figure [6} taken from Figu re El 
in Appendix 5 from |Epinat, Amram &; Marcelin|2008 ), they 
show some common features when compared to the rotation 
curve from GHaFaS data: a dip of the curve in the central 
part of the galaxy, then a steep rise up to 30 arcsec from 
the centre (i.e., up to 3.04 kpc from the centre), and a less 
steep slope in the last part of the curve. Some of these fea- 
tures found in rotation curves from Ha observations can be 
a consequence of the reduction process (e.g., wrongly fixed 
parameters when deriving the rotation curve, such as the 
systemic velocity), or can be caused by a physical feature 
(e.g., the presence of a bar) or even due to the observations 
(patchiness in the line emission). 

4-2.2 Velocity models and residual maps 

The presence of non-circular motions in spiral galaxies can 
have a considerable impact on the shape of the rotation 
curve, and thus on the amount of dark matter inferred. On 
the other hand, the non-circular motions carry important 
information about the components that cause them. For in- 
stance, the strongly non-axisymmetric gravitational poten- 
tial of the bar can cause local deviations on the velocity field 
(streaming motions on the gas flow of galaxies have been dis- 
cussed by, e.g., 



Duval & Athanassoula 1983 under the beam 



scheme; by Lindblad, Lindblad & Athanassoula 1996 using 
the FS2 code presented in|van Albada &; Roberts|1981| |van| 
Albada, van Leer & Roberts 1982 and van Albada 1985; or 
by jSalo et al.|1999| using the sticky-particle method), while 
strong spiral structure can provoke shocks and streaming 
motions in the gas (e.g., Bosma|1978| Visser|1 980; Marcelin, 



|Boulesteix 8z Georgelin 1985; R ots et al.| 1990; Kn apen et al 
|l993| and |Knapen, Shlosman fe Peletier|2000| ) . 

To study the non-circular motions in NGC 864, we use 
the rotation curve derived above from the GHaFaS data to 
construct a model, assuming kinematic symmetry, and us- 
ing the GIPSY task VELFI. This procedure creates a model 
velocity filed plotting the information from the provided ro- 
tation curve into the same dimensions as the original map 
(Figure [7]). A residual velocity map (Figure [sj is created by 
subtracting the model from the original velocity map. This 
residual map can be interpreted as a view of the non-circular 
motions. 

These figures show that the region with the highest non- 
circular motions, of around 50 km/s in absolute value (about 
one third of the rotational velocity), is delineated by the bar. 
As explained before, it is easier to measure and highlight 
these non-circular motions in the smoothed map. 

4-2.3 Position- velocity diagrams 

One can study the non-circular motions by examining the 
residual velocity field, as we saw in the previous section. 
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Figure 8. Residual velocity maps obtained from the low-resolution map {left) and from the high-resolution map {right). Bottom: the 
bar zone for both maps is highlighted. 



However, as each point on the rotation curve is calculated 
by azimuthally averaging the complete velocity information 
(including non-circular motions) within a certain range, the 
rotation curve and the resulting model velocity field will be 
influenced to some degree by non-circular motions. Conse- 
quently, we used position- velocity diagrams (PV diagrams), 
because they reproduce all the observed information, even 
though they are not always easy to interpret. 

We used the SLICE task in GIPSY to extract a plane 
from a data set at the angles of the kinematic major axis 
(20°), the kinematic minor axis (110°) and along an inter- 
mediate angle of 70° (wedges of ±5° were used in all cases). 
The data used are the low-resolution GHaFaS cube, and the 
resulting PV diagrams are presented in Figure [9] 

A PV diagram of a galaxy along the kinematic minor 
axis, in this case 110°, is completely flat in the absence of 
non-circular motions, whereas in the same ideal case a PV 
diagram along the kinematic major axis shows the rotation. 
At intermediate angles, such as the 70° slice we show for 
NGC 864, an intermediate behaviour is expected. In the case 
we present here, the PV diagrams show clearly the presence 
of non-circular motions, and indicate their location and ap- 
proximate amplitude. We will describe in the next section 
that they correspond to motions induced by the bar and 
spiral arms. The centre shows a very broad extent of veloc- 



ities in Ha emission, suggesting a strong velocity gradient 
here. Additional contributions from strong star- forming ac- 
tivity in the nucleus, and/or incorrect continuum subtrac- 
tion there, cannot be completely excluded. We will ignore 
the centre in our further analysis. The minor axis PV dia- 
gram shows a feature at 2^15 m 29 s and ^1660 km/s. It arises 
from the [NIIJ-6548 line leakage that didn't affect the com- 
putation of the moment maps, but that shows up in a PV 
diagram. 



4.3 Star Formation Rate 

The star formation rate (SFR) was determined following 
Kennicutt et al. (2009) with the expression: 



SFR(M yr" 1 ) = 5.5 x 10" 42 L(Ha) (1) 
where L{Ha) is the luminosity, calculated as 



L(Ha)[erg/s] = 4ttL> 2 (3.086 x lO 24 ) 2 ^, 



(2) 



with D the distance to the galaxy in Mpc (Table 1) 
and F^ a the flux corrected for Galactic and internal ab- 
sorption. In this case, the Galactic absorption was taken 
from NED, with a value of A(R) = 0.130 mag (calculated 
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Figure 9. PV diagrams for NGC 864 from the low-resolution 
GHaFaS velocity data, a) PV along the kinematic major axis 
(PA=20°, counter-clockwise from N); b) PV along an interme- 
diate angle (70°) and c) PV along the kinematic minor axis 
(PA=110°). Contour levels are at 0.5, 1.5, 4.5, 13.5 and 25<j. 
Overlaid (grey dots) is the velocity profile for the whole disk at 
the corresponding angle and the systemic velocity (dash-dot-dash 
line) corresponding to a velocity of 1562 km/s. 



Table 2. L(Ha) and SFR measured from the AC AM Ra image 
for the whole galaxy, the eastern and western arms, for the bar 
region, and for the brightest regions named with A, B, C and D 
identified in panel a) of Figure [2] (the centre labelled with A; B 
in the south, at the end of the western arm, and C and D at the 
eastern and western ends of the bar, respectively). 



Region 



L(Ra) (erg s" 1 ) SFR (M© yr" 1 ) 



Whole galaxy 


(4.0 ± 


1.6) 


10 41 


2.19 


± 


0.88 


Eastern arm 


(1.5 ± 


0.6) 


10 41 


0.85 


± 


0.34 


Western arm 


(1.9 ± 


0.8) 


10 41 


1.04 


± 


0.41 


Bar 


(5.3 ± 


2.1) 


10 39 


0.03 


± 


0.01 


A (centre) 


(3.7 ± 


1.5) 


10 40 


0.20 


± 


0.08 


B (End of western arm) 


(2.0 ± 


0.8) 


•10 40 


0.11 


± 


0.04 


C (start of eastern arm) 


(1.0 ± 


0.4) 


.1 40 


0.06 


± 


0.02 


D (start of western arm) 


(3.8 ± 


1.5) 


.1 39 


0.02 


± 


0.01 



based on the dust maps frompchlegel, Finkbeine r &; Davis| 
1998). To correct for internal absorption, there are several 
proposals in the literature (see Sanchez- Gall ego et al.||2012 
for an overview). We have adopted the constant value of 
A(Ha) = 1.1 mag ( |Kennicutt Kent||1983| ). 

We have measured the SFR in the brightest regions 
(marked with letters in panel a) of Figure [2]), and also for 
the whole galaxy. The final values are presented in Table 
[2] Not all the Ha emission in the centre is necessarily due 
to massive star formation, with contributions possible from 
moderate non-stellar activiy (even though NGC 864 appears 
not to have an AGN activity, see |Ho, FiTi ppenko & Sargent 
|1997| ) or from shocks. This value should thus be considered 
as an upper limit, of around one tenth of the total SFR in the 
galaxy. The SFR for the bar region is low (0.03 M yr _1 ), 
which can be explained by a dominant old population there. 
There is some massive star formation from the regions at 
the end of the bar, but almost all massive star formation in 
this galaxy occurs in the two main spiral arms. 

There are several factors to take into account to es- 
timate the uncertainties in the measurements of the SFR. 
First of all, the adopted value for the galaxy distance taken 



from NED refers to Tully et al. (2009), which used the Tully- 



Fischer relationship to compute the galaxy distance. There- 
fore, a typical 20% uncertainty is estimated for the distance 
value. The basic reduction processes, i.e. flat-fielding cor- 
rection, can contribute with around 2% of uncertainty. Sec- 
ondly, the zero-point calibration carries an uncertainty of 
3%, in agreement with the value of 2% typical for photo- 
metric nights. The uncertainty in the flux measurement is 
around 1%. However, apart from the uncertainty due to the 
distance, the highest uncertainty comes from the continuum- 
subtraction, estimated to be around 11%. The scaling factor 
is somehow a subjective value, and to have an idea of the 
uncertainty, we have measured the scaling factors from two 
images: one slightly oversubtracted and another slightly un- 
dersubtracted, leading to an uncertainty of the 11%. 

Altogether, taking into account the uncertainties re- 
lated to the quality of the image (without considering the 
distance) we estimate an uncertainty of 17%. The estimated 
uncertainty is comparable to previous studies using the same 



method: James et al. (2004) who found a value of around 



14%, Kennicutt et al.l (12008) who found a value of 16%, and 



Sanchez- Gallego et al.|(2012) who found a value of 18 % 
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5 DISCUSSION: INFLUENCE OF THE BAR 
ON THE KINEMATICS 

The analysis of non-circular motions is of great interest to 
the understanding of galaxy structure and kinematics. The 
rotation curves (Figure |6|, PV diagrams (Figure [9]), and the 
residual maps (Figure pT provide crucial information about 
the kinematical processes that take place in the galaxy. A 
prominent bar, such as the one in NGC 864, is often the 
most important factor in this context, and indeed the pres- 
ence of the non-axisymmetric gravitational potential of the 
bar can be recognised in the deviations from circular mo- 
tion across the galaxy. In the residual maps (Figure |8j), the 
most prominent non-circular motions are found along the 
bar and at the beginning of the eastern arm. This confirms 
that the bar has a significant influence on the kinematics of 
the galaxy, causing the velocities to deviate from the circu- 
lar, rotational, motion. The kinematics along the kinematic 
minor axis may further quantify the influence of the poten- 
tial of the bar, due to the fact that the PA of the bar is 
very similar to the PA of the kinematic minor axis. Fig- 
ure [To] thus presents the GHaFaS velocity profile along the 
kinematic minor axis (top panel), and also the shape of the 
isovelocity contour at the systemic velocity (lower panel). 

The profile shows that there are significant changes in 
the velocity along the kinematic minor axis, of up to 50 
km/s. In the profile, it is noticeable that most of the veloc- 
ities along the kinematic minor axis are below the systemic 
velocity. Various peaks and troughs are found along the pro- 
file. In order to understand why, and avoiding at the same 
time the problem of the patchiness of the Ka emission, we 
repeat this exercise with the SAURON H/3 and stellar veloc- 
ity maps. Here, the peaks and troughs can easily be distin- 
guished, and the deviations from the systemic velocity also 
range up to 50 km/s. 

One should notice the differences between these 
SAURON maps and the one from GHaFaS. Firstly, the FOV 
of SAURON (33" x 41") is much smaller than the one of 
GHaFaS (202" x 202"), so the profiles along the kinematic 
minor axis in this Figure are only the central part of the 
profile obtained with GHaFaS data. Secondly, the patchi- 
ness implies a lack of points in some parts of the GHaFaS 
profile. In Figure [TT] we present the bar region of NGC 864 
as seen in the low-resolution velocity map from GHaFaS 
data, and overlying, the intensity contours from SAURON 
H/3 observations and the contours of the derived systemic 
velocity from SAURON H/3 and stellar velocity maps. 

In the profile from the GHaFaS velocity map, there are 
lower velocities in the eastern part of the bar, and higher ve- 
locities in the western part. PV diagrams (Figure [9]) confirm 
this. The PV diagram along the kinematic major axis shows 
the rotation, whereas the PV diagram along the kinematic 
minor axis confirms the deviations along the bar with the 
peaks mentioned above. The correspondence between peaks 
can indicate bar streaming motions, although these results 
are not nearly as conclusive as in, for example, the much 
clearer case of M100, as described in, e.g., |Knapen, Shlos-| 



1600 



man & Peletier 



(2000). 



Taking everything into account, the profiles and PV di- 
agrams along the kinematic minor axis show that there are 
deviations from the circular rotation velocity, confirming the 
presence of non-circular motions along the bar. These are 
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Figure 10. Top: the profile along the kinematic minor axis of 
NGC 864 in black dots; overlaid in a blue dash-dot-dash line the 
systemic velocity of 1562 km/s; overlaid in red crosses and green 
squares the 1562 km/s profiles at kinematic minor axis from the 
SAURON ¥1(3 velocity map and stellar velocity map. Bottom: 
velocity map obtained with GIPSY. The red dash-dot-dash line 
shows the location of the kinematic minor axis, whereas the black 
contours trace a velocity of 1562 km/s, corresponding to the sys- 
temic velocity. In grey, the location of the SAURON image to 
compare the different scales. 



probably caused by the non-axisymmetrical potential cre- 
ated by the bar, as widely discussed in the literature (see 
references in Section 4.2.2). 



6 CONCLUSIONS 

We have presented a kinematic study of the spiral galaxy 
NGC 864. Data sets obtained with three different WHT in- 
struments have been studied, resulting in rotation curves 
and intensity, velocity, residual and gradient maps. From 
this, we reach the following conclusions: 

(i) A technique for flux calibration of GHaFaS data has 
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Figure 11. Bar region of the low-resolution velocity map of NGC 
864 from the GHaiFaS data. Overlaid are contours from SAURON 
data: the dashed thick line corresponds to a velocity of 1562 km 
s _1 in the H/3 velocity maps, whereas the dotted line corresponds 
to a velocity of 1562 km s — 1 in the stellar velocity map. Overlaid 
black thin line corresponds to 1562 km/s in GHaPaS data. 



been developed. This will be used for other S 4 G galaxies 
observed with GHaFaS. 

(ii) Our bulge/disk/bar 2D decomposition of the S 4 G 
3.6/ira image of NGC 864 using budda yields a flat-profile 



bar, which is typical of strong bars ( |Athanassoula &; Misirio- 
tis|2002 ) and of two-armed grand design spirals ( |Elmegreen 



fc Elmegreen|1985|). Flat-profile bars are linked to a higher 



degree of exchange of angular momentum with the halo than 
exponential bars ( Athanassoula|2003 ) . Such bars can result 
from sharply rising rotation curves, as observed in NGC 864, 
and are most likely connected to the observed spatial offset 
between the ends of the bars and the onset of the spiral 
arms. 

(iii) The rotation curves obtained from the GHaFaS ve- 
locity maps are similar to those obtained by |Epinat, Amram] 
|fe Marcelin] ( |2008| ), both using Ha observations. However, 
the curves do not cover all the disk and the intrinsic patch- 
iness linked to the line emission prejudices the analysis and 
the derived results. Due to the high angular and spectral 
resolution, our new rotation curve shows considerable detail 
in the central 2 kpc radius. 

(iv) We have found non-circular motions along the bar in 
the residual maps, in the PV diagram along the kinematic 
minor axis, and in the velocity profile along the kinematic 
minor axis. These velocity patterns are typical of a barred 
galaxy behaviour, with streaming motions along the bar. 

This paper is the first in a series presenting the results 
of the programme to observe a sample of some 40 spiral 
galaxies of all morphological types from the S 4 G sample with 
the GHaFaS FP instrument. 
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